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Kepler’s Supernova: An Overluminous Type la Event Interacting 
with a Massive Circumstellar Medium at a Very Late Phase 
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ABSTRACT 

We have analyzed XMM-Newton, Chandra, and Suzaku observations of Ke¬ 
pler’s supernova remnant (SNR) to investigate the properties of both the SN 
ejecta and the circumstellar medium (CSM). For comparison, we have also ana¬ 
lyzed two similarly-aged, ejecta-dominated SNRs: Tycho’s SNR, thought to be 
the remnant of a typical Type la SN, and SNR 0509-67.5 in the Large Magellanic 
Cloud, thought to be the remnant of an overluminous (SN1991T-hke) Type la SN. 

By simply comparing the X-ray spectra, we hnd that line intensity ratios of iron- 
group elements (IGE) to intermediate-mass elements (IMF) for Kepler’s SNR 
and SNR 0509-67.5 are much higher than those for Tycho’s SNR. We therefore 
argue that Kepler is the product of an overlumious Type la SN. This inference 
is supported by our spectral modeling, which reveals the IGE and IME masses 
respectively to be 0.95to;37MQ and 0.12to;o5MQ (Kepler’s SNR), 0.75to;®5MQ 
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and 0.341 °;°®Mq (SNR 0509-67.5), and 0.35l[(;^^ Mq and O.7Ol)J;^^M0 (Tycho’s 
SNR). We find that the CSM component in Kepler’s SNR consists of tenuons 
diffnse gas (~O.3M0) present thronghont the entire remnant, plus dense knots 
(~O.O35M0). Both of these components have an elevated N abundance (N/H~4 
times the solar value), suggesting that they originate from CNO-processed mate¬ 
rial from the progenitor system. The mass of the diffuse CSM allows us to infer 
the pre-SN mass-loss rate of the system to be ~ 1.5x10 ®(n„/10kms ^)M0yr 
in general agreement with results from recent hydrodynamical simulations. The 
dense knots have slow optical proper motions as well as relatively small X-ray- 
measured ionization timescales, which indicates that they were located a few 
pc away from the progenitor system and were only recently heated by forward 
shocks. Therefore, we argue that Kepler’s SN was an overluminous (91T-like) 
event that started to interact with massive CSM a few hundred years after the 
explosion. This supports the possible link between 91T-hke SNe and the so-called 
“la-CSM” SNe — a rare class of SNe la associated with massive CSM. The link 
implies that ~10% of SNe la are associated with massive CSM which most likely 
originates from a companion star in a single degenerate progenitor system. 

Subject headings: ISM: individual objects (Kepler’s SNR, Tycho’s SNR, SNR 0509- 
67.5) — ISM: supernova remnants — circumstellar medium — supernovae: gen¬ 
eral — X-rays: general 


1. Introduction 


Type la supernovae (SNe la) are important astrophysical objects, as they play a central 
role in measuring the accelerative expansion of the universe and are primary source of iron- 
group elements (ICE) in the universe. However, the progenitor syste m(s) and detailed 


explosion mechanism(s) have been a matter of considerable debate (e.g., iMaoz et al.ll20l4 
for a recent review). It has been believed that SNe la are runaway thermonuclear explosions 
of a C/0 white dwarf (WD), as it approaches the Chandrasekhar mass (~1.4M0). There 
are mainly two long-standing competing scenarios for how a WD increases its mass: (1) the 
single-degenerate model, in which a WD accretes rn aterial from a non-degenerate companion 
star in a close binary system (I Whelan fc Ibenlll97311: and (2) the dou ble-degenerate model, in 


which two C/0 WDs merge fllben fc Tutukovi 11984 IWebbinklll984l) . Therefore, identifying 


the progenitor stars (as a means of discriminating between the two scenarios) is an area of 
considerable activity in SN la studies. 


There are many observational approaches to identifying progenitors, including studies 
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of potential progenitor populations flHachisu fc Katoll200lh . sea rches for companions to the 


2014; 

McCullv et ah 

2014) 

and surviving companions in nearby supernova remnants (SNRs) 

(e.g., 

Schaefer &: Pagnotta 

2012; 

Bedin et ah 

2014 

), and evidence of interactions between the 


SN ejecta and the companion star fe.g.. lMaeda et al.ll2014j) . In addition, recent detection of 
a circumstellar medium in some SNe la (CSM: environmental gas modihed by the mass loss 
from the progenitor systems) have opened a new window to infer progenitor systems. 

The presence of a H-rich CSM itself supports the SD scenario, since such a CSM is not 
generally expected around two WDs in the DD scenario. Analysis of the CSM signature (both 
absorption and emission) offers further insight into the progenitor system. For example, time- 
variable Na I D absorption features, which were found in the normal T ype la SN 2006X and 
a few others flPatat et al.ll2007l : iBlondin et al.ll2009l: ISimon et al.ll2009l) . could be interpreted 
to arise frorn expa nding s hells (or knots) due to s uccessive nova eruptions of the progenitor 
fjPatat et al.l 1200711 . Later, ISternberg et al.l (120111) found that about half of all sampled SNe 
la (12/22) in their sample show blueshifted Na I D absorption lines, while only a quarter 
of them (5/22) show redshifted absorption, and the remaining quarter (5/22) show single or 
symmetric absorption prohles. They suggested that the biased detection of the blueshifted 
systems means a large fraction of SNe la (^25%) are associated with progenitors’ gas outflows 
(otherwise, we would see the same fractions of blueshifted and redshifted prohles). From the 
absorption column densiti es ^ 10^^cm ~^), the masses of the CSM shells have been 


estimated to be <0.01 Mq flPatat et al.ll2007l) . which are consistent with successive recurrent 
nova eruptions that swept up the stellar wind of the companion. Therefore, these CSM 
properties can be best explained by a symbiotic recurrent nova progenitor, similar to RS 
Ophiuchi. 

Evidence for much more massive CSM (~O.O1-5M0) has been found in several SNe 
la (SN 200 2ic, SN 2005gj, SN 1997cy, SN 1999E, S N 2008J, PTFllkx, and some other 


candidates: iHamuv et al.l 120031: ISilverman et al.lI2013L references therein). These SNe are 


called “SNe la-CSM” which are characterized by a narrow Ha line on top of an overluminous 
(SN 1991T-hke or 91T-like) la te mplate in the optic al spectra. SNe la-CSM comprise only a 
small fraction of SNe la: 0.1-1 % flDildav et al.ll2012[) . albeit a true fraction may be somewhat 
larger, given a possibly signihcant amount of SNe la-CSM that are incorrectly classihed as 
Type Iln because of their spectral resemblance ( Leloudas et al.lbolO ). While there are several 
proposed progenitor systems for la-CSM objects, none of them is yet conclusive. 

Kepler’s SNR, the remains of SN 1604 A.D., is one of a handful of SNRs which are associ¬ 
ated with historical Galactic SNe. Located at only 3-7 kpc flKerzendorf et all 2014 . references 
therein), this remnant provides a precious opportunity to study details of SN la explosions 
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and the progenitor system. It belongs to a rare class of Type la SNRs flKinugasa fc Tsnnemi 
19991: iRevnolds et al.l 120071) that exhibit emission from a CSM, evidenced by the extremely 
high ar nbient density (^100 cm"^) at >500 PC o nt of the Galactic plane and a N overabnn- 
dance flBlair et alJll99ll: iGerardv fc Fesenll200ll) . Optical monitoring of 50 long-lived GSM 
knots revealed that they are moving fairly slowly, wit h a space velocity of r'U 280 km s ^ at a 
distance of 4.5 kpc (IBandiera fc van den BerghI 1199111 : reminiscent of quasi-st at ionary floc- 
cnli in Gassiopeia A. The trace-back time was measured to be 3.2 ± 1.2 x lO^yrs, which 
should be taken as a lower limit, given possible acceleration by the SN forward shock. 


Douvion et al.l 0200ll) found infrared (IR) emission from warm dust. Since the IR inten¬ 


sity map is similar to the distribution of the dense optical CSM knots, the au thors argued 
that the IR emission arises from dust in the shocked CSM. iBlair et al.l (120071) further sup¬ 
ported this picture based on Spitzer Space Telescope observations, and estimated the total 
shocked CSM mass to be ~1 Mq from their relatively cursory analysis of X-ray data. In 
addition to the CSM knots, t he presence of a diff use CSM has been inferred using fol¬ 
lowup Spiteer observations bv IWilliams et al.l (120121) . Their work indicates that the CSM 
dust consists of two components: lower-temperature (~80-100 K) dust behind fast shocks 
(>1000kms“^) and penetrating into the tenuous ambient medium produces the majority 
of the IR emission; higher-temperature (>150 K) dust behind slower shocks (a few hundred 
kms“^) penetrating into high density material {n ~ 50 — 250 cm“^) coincides with the optical 
dense k notslll The presence o f a CSM is also supported from recent hydrodynamic simu¬ 


lation s ( Velazquez et al.l l2006l: IChiotellis et al.l l2012l: iToledo-Rov et al.l l2014l: iBurkev et ah 
20131) . where a high mass-loss rate of M ~ 10“® Mq yr“^, suggestive of dense winds from an 


asymptotic giant branch (AGB) star, can reasonably explain the observed expansion index 
and the north-south asymmetric morphology of Kepler’s SNR. 

X-rays arise predominantly from SN ejecta, which are not detected at either optical or 
IR wavelengths. The X-ray spectrum exhibits strong Fe L lines and weak O K lines. A 
multi-component non-equilibr ium ionization modeling sug gests a Type la SN rather than a 
core-collapse SN as its origin (iKinugasa fc T.snnemilll999L references therein for a history of 
X-ray observations). From deep Chandra observations, IRevnolds et al.l (120071) conhrmed the 
high Fe /0 ratio and fo u nd no evidence for a neutron star, further supporting the Type la 
origin. iPatnaude et al.l (120121) claimed that Kepler’s SN was likely to be a 91T-like event 


that produced ~1 Mq of ®®Ni, by comparing the observed X-ray spectrum with simulated 
spectra bas e d on their hydrodynamic modeling at 400 yrs after the explosion. More recently. 
Park et al.l (120131) performed a deep Suzaku observation, and found a high Mn/Cr line in- 


^We note that dust is mainly heated by collisions in shocked gases, and that the efficiency of collisional 
heating increases with the density of the shocked gas. 
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tensity ratio of ~0.6, leading them to snggest a snpersolar metallicity (~3 times the solar 
valne) of the progenitor star. 

Here, we investigate X-ray spectra of Kepler’s SNR and two other yonng la SNRs, 
Tycho’s SNR and SNR 0509-67.5, for comparison, by making the best nse of available data 
from currently operational X-ray observatories. In the following section 2, we summarize the 
observations used. In sections 3 and 4, we perform X-ray spectral analyses, focusing mainly 
on the SN ejecta and the CSM, respectively. Finally, we discuss the observational results 
and lead to conclusions. 


2. Observations and Data Reduction 


We use archival data obtained by XMM-Newton, Chandra, and Suzaku, as summarized 
in Tabled! While our object of interest is Kepler’s SNR, we also analyze two similarly-aged la 
remnants, Tycho’s SNR and SNR 0509-67.5 in the Large Magellanic Cloud (LMC). Detailed 
classifications have already been determined for these tw o remnants, based o n their light-echo 
spectra: Tycho is thought to be a prototypical Typ e la flKrause et al.ll2008[) . SNR 0509-67.5 
is thought to be “SN 1991T-hke” flRest et al.l 120051) . Thus these two remnants can serve as 
ideal reference targets to study the amount of IGE in Kepler’s SN. 


We reprocessed the raw data, using version 13.5.0 of XMM-Newton Science Analysis 
Software, version 4.6 of Chandra Interactive Analysis of Observations, and version 22 of the 
Suzaku software, together with the latest ver sions of the calibratio n files at the analysis phase. 
For the X-ray Imaging Spectrometer (XIS: iKovama et ahl 120071) onboard Suzaku, we focus 
on front-illuminated CCDs (XISO and XIS3) which have slightly better spectral resolution 
than the back-illuminated CCD (XlSlj. We generated respons e files for the XMM-Newton 
Reflection Grating Spectrometer IRGS: Iden Herder et al.ll200ll) . taking account of energy- 
dependent spati al structures in the source, as we did in our previous analyses of the Galactic 
SNR Puppis A ( Katsuda et ah 2012 . 2013h . The hrst-order RGS spectral resolution, which 
is determined by the size of the remnants in our case, is roughly E/AE ~25 and 180 for 
Kepler’s SNR and SNR 0509-67.5, respectively. This resolution is higher than that of the 
XIS of E/AE ~15. 


Backgrounds (BGs) were taken from source-free regions in the same fields of view. The 
exception is the RGS spectrum for Kepler’s SNR, for which the spatial extent of the source- 
free areas is limited; hence we used blank-sky data (Obs. ID; 0051940401) taken immediately 
after the observation of Kepler’s SNR. Each spectrum was grouped into bins with at least 
20 counts in order to allow us to perform a f^st. 
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3. Ejecta Masses from Integrated X-Ray Spectra 


3.1. Glancing at Suzaku'’s XIS Spectra of the Three SNRs 


We first examine the wide-band X-ray spectra for the three SNRs. Figure [T] shows 
integrated Suzaku XIS spectra of the SNRs, where the intensities are normalized to equalize 
the Si Hea line. A number of emission lines are apparent, mostly arising from the shocked 
SN ejecta. Obviously, the strengths of Fe lines (both L- and K-shell transitions) are quite 
different among the three sources, whereas K-shell lines from intermediate mass elements 
(IMF: Si, S, Ar, and Ca) are similar to one another, in the sense that shock-heated IGE- 
to-IME mass ratios for Kepler’s SNR and SNR 0509-67.5 are signihcantly larger than that 
for Tycho’s SNR. It is reasonable to consider that the amount of unshocked cold ejecta in 


evolutionary states ( 

Dickel et al. 

1988; 

Revnoso et al. 

1997; 

Katsuda et al. 

2008; 

^ -J ^ / 

Vink 

2008; 

Katsuda et al. 

2010 j. and the similar reverse shock positions f~70% of the forward shock 


radius) as we describe below. Therefore, Kepler’s SN must have produced a relatively large 
mass of IGE, similar to SNR 0509-67.5. When combined with information from light-echo 
spectra, i.e., overluminous for SNR 0509-67.5 and normal for Tycho’s SNR, we argue that 
Kepler’s SN was likely an overluminous la event. 


This argument will be substantiated by our spectral htting in the following section. 
There are reasons why spectral htting is required to estimate the ejecta masses quantitatively 
and correctly. For example, while SNR 0509-67.5 has the strongest apparent Fe L lines below 
IkeV in the raw XIS spectra in Fig. [H these lines are considerably weaker than those of 
Kepler’s SNR, if we correct for the interstellar absorption. Also notable is the slight shift of 
the line center energies of SNR 0509-67.5 with respect to those of the two Galactic SNRs. 
This shift is due not to artihcial/instrumental effects, but to a mixture of bulk motions and 
ionization effects, as we will show below. Since the ionization state is sensitive to X-ray line 
intensities, we need to care about it in estimating ejecta masses. All of these effects can be 
accounted for through spectral htting. 


3.2. Quantitative Mass Estimates Based on Spectral Fitting 

In our spectral htting, we use XMM-Newtorks RGS below 2 keV and Suzaku^s XIS above 
2keV for Kepler’s SNR and SNR 0509-67.5. On the other hand, for Tycho’s SNR, we use 
only data from Suzaku^s XIS, since the “slitless” RGS cannot provide a high-resolution spec¬ 
trum due to the relatively large angular size of the remnant (diameter of ~8.5'). Spatially- 
integrated X-ray spectra for the three SNRs are shown in Fig. [2l In these plots, only the 
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first-order RGS spectra are shown, but we use both hrst and second orders when modeling 
the data. The data points in blue for Kepler’s SNR is a local-BG subtracted GSM spectrum 
extracted from the summed-up regions of knots 1 and 2 (see, Section 4), by using Chandra’’s 
Advanced GGD Imaging Spectrometer (AGIS: iGarmire et ahN2003l) . We fit this spectrum 
with a GSM-only model having all normalizations of ejecta components fixed to zero, so that 
we can separate out the GSM emission as cleanly as possible. Note that, for the spectrum 
of Tycho’s SNR, we exclude data below 0.6 keV owing to poor signal-to-noise ratios due 
to a relatively heav y interstellar absorp tion together with the low energy tail emission from 
higher energy lines flKovama et ahll2007l) . and in the 1.55keV-1.75keV band due to a Si-edge 
calibration concern with the xid§. 

We fit these spatially-integrated spectra with an absorbed, vpshock -f three vneis + 
power-law -f- several Gaussian components model in the XSPEG package fjArnaudlll996l) . 
Here, the vpshock and the power-law components are supposed to be the thermal and 
synchrotron radiation, respectively, from the swept-up medium. The three vnei components 
represent emission from SN ejecta. At least three vnei components are required from a 
statistical point of view (F-test probability of >99%); we chose to not add more in order to 
keep our model as simple as possible. The several Gaussian components represent lines of 
Fe L and/or O K, Fe L and/or Ne K, Gr K, and Mn K not well represented in the broad 
band models. 

While details of our parameter treatment can be found in Table [2], we elaborate some 
of them below. We £x the abundance of either Si/H or Fe/H to be 10® times the solar value, 
assuming that the three vnei components are nearly pure metal plasmas. Most of the emis¬ 
sion from IMF is reproduced by two vnei components, for which the two-components nature 
may reflect a temperature/density gradient within the shocked ejecta; e.g., a hydrodynamical 
model of Type la SNR evolution in a uniform ambient density predicts a strong temperature 
decrease (as well as de nsity increase) from the reverse shock to the contact discontinuity 
fjPwarkadas et al.lll998l) . For these components, we tie the metal abundances to each other 
except for Fe (and Ni). The remaining vnei component is introduced to reproduce most of 
the Fe K emission. We as sume the Ni/Fe abun dance ratio to be that expected in the classical 
deflagration model (W7: iNomoto et al.l 119841) : the assumed value of Ni/Fe abundance does 
not affect our spectral fitting, since Ni lines are not clearly detected in the spectra. In this 
Fe-rich component, we allow abundances of Ar and Ga to vary freely; otherwise we see sig¬ 
nificant residuals at low-energy sides of the K lines of these elements, especially for Tycho’s 
SNR. This would indicate that some Ar and Ga may be mixed into IGF-rich layers, although 


^http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/sical.html 


















we cannot rnle ont other possibilities, e.g., the presence of strongly redshifted components 
for these elements. In any case, the ejecta masses will not change dramatically, hence we 
here adopt the simplest solntion in terms of the spectral modeling. 


We fix the absorbing colnmn densities at iVn = 6.4 x 10^^ cm~ lx 10^^ cm and 
6 X 10^° cm“^ adopting elemental abundances given bv I Wilms et ahl (120001) for Kepler’s SNR, 
Tycho’s SNR, and SNR 0509-67.5, respectively. The values for the two Galactic SNRs are 
determined by fitting synchrotron-dominated spectra from the outermost thin rims obtained 
by the Chandra's ACIS, while that for SNR 0509-67.5 is taken from the literature (e.g., 
Warren et ahl 12004 IWilliams et al.ll201l[) . Free parameters are the electron temperature. 


kTe, the ionization timescale, Ugi, the normalization dehned as J UenudV where Ue and 
riH are electron and hydrogen densities and V is the emitting volume, and the redshift of 
each component, some of which are linked among different components in order to better 
constrain the parameters. Photon indices are allowed to vary freely. In addition, we multiply 
the gsmooth model to each thermal component in order to take account of line-broadening 
effects, for which we assume the widths to be proportional to the line energies. 


As shown in Fig. [21 the best-fit models in green represent the data well. Individual 
components are shown in blue (swept-up medium), red (ejectal = the lower-temperature 
IME-dominated component), orange (ejecta2 = the higher temperature IME-dominated 
component), light blue (ejectaS = the Fe-rich component), gray (power-law), and dotted 
black (Gaussians). Note for Kepler’s SNR that the swept-up components seen for both 
XMM-Newton's RGS and Chandra's AGIS have the same intrinsic spectral shapes, but look 
different due to distinct spectral responses. The lower panels in Fig. |2]show residuals; the 
red points represent the RGS second order data that are not shown in the upper panels. The 
reduced-y^ values summarized in Table [2] are 1.38 (SNR 0509-67.5), 1.45 (Kepler’s SNR), 
and 2.91 (Tycho’s SNR). The relatively large value for Tycho’s SNR is not surprising, given 
that it is the brightest source of the three SNRs and the larger spectral bin widths (3.65 eV 
at the minimum using the XIS) than those of the other two sources (<1 eV using RGS), mak¬ 
ing systematic uncertainties relatively serious for this SNR. We do not, however, introduce 
systematic uncertainties in our spectral fitting, since they do not affect our main conclusions. 


We check the robustness of our spectral analysis by comparing the best-fit parameters in 
Table [2] with previous measurements. First of all, the N abundance of the swept-up (vpshock) 
component in Kep ler’s SNR is quantitatively consistent with previous optical measurements 
(IBlair et al.lll9911 1. showing that this component is the GNO-processed stellar wind material, 
i.e., the GSM, as we expected. Also, the lines in the GSM component (mainly N Lya, O Hea, 
and O Lya) are signihcantly narrower than those of ejecta components, assuring robustness 
of the spectral separation between the GSM and the ejecta. Unfortunately, the data do 
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not allow us to measure N abundances or line widths of the CSM components for the other 
two SNRs because of the relatively poor spectral resolution (Tycho’s SNR) and statistics 
(SNR 0509-67.5). A s for the photon indic e s, the best-ht values agree with previous X-ray 
measurements (e.g., ITamagawa et al.l 120091: ICassam-Chenai et al.l 120071 : IWarren et al.l I2004J: 
Kosenko et al.l 120081) . In addition, by converting the line broadening of the ejecta compo¬ 
nent to bulk velocities, we obtain = 2740^6° km s“^, 2970tio kms“^, ^^^1 3570^230 km s“^ 
for Kepler’s SNR, Tycho’s SNR, and SNR 0509-67.5, respectively, where the velocities in¬ 
clude effects of both thermal broadening and shell expansion whose contributions are hard 
to distingui sh from our data. T hese are somewhat s maller than previous estimates for Ty¬ 
cho’s SNR ( Havato et ah 2010 ) and SNR 0509-67.5 ( Kosenko et ah 2008 ). which might be 
explained by multiple ionization effects in our analysis. 

We next look into details of ejecta components. Figure [3] illustrates unabsorbed spectra 
of the total ejecta components with intensities normalized at Si K-shell lines, in order to 
clearly reveal the differences among the three SNRs. We can readily see that the Fe K and 
L lines for Kepler’s SNR and SNR 0509-67.5 are much stronger than those for Tycho’s SNR. 
This is also evident from the spectral ratios in the lower panel. We then quantitatively 
estimate ejecta masses, using the best-ht parameters listed in Table [H As the hrst step, 
we calculate Ue/riH for each component, t aking account of the elemental abundances and 
the ionization balance (in the SPEX code: iKaastra. et al.l 119961) at the fcTg and nj, values 
measured. With this information, we can calculate J which can be converted to 

j n\dV by multiplying each elemental abundance. 

To deduce densities (and masses) from / n^dV, we need to know emitting volumes 
of individual components. To this end, we generate deprojected radial prohles from high- 
angular resolution data obtained with the Chandra’s ACTS. We hrst extract energy-band 
images for 0.7-0.85keV (Fe L), 0.9-1.2keV (Fe L), 1.65-2.0keV (Si K), 2.25-2.6keV (S K), 
2.95-3.25keV (Ar K), 3.65-4.1 keV (Ca K), 4.0-6.0keV (continuum), and 6.1-6.8keV (Fe 
K). For each band image, we subtract continuum emission extrapolated according to the 
power-law slopes given in Table [2l We then generate continuum-subtracted radial prohles, 
by focusing on a smooth limb-brightened region of each SNR, i.e., north, west, and north¬ 
east for Kepler’s SNR, Tycho’s SNR, and SNR 0509-67.5, respectively, as shown in Fig. H] 
upper panels. The radial prohles produced in this way for four selected energy bands are 
displayed in the middle panels of Fig. 01 We see a general trend, with Fe K, Fe L, and Si 
K pea king at increasing radius, co nsistent with a previous XMM-Newton study for Kepler’s 
SNR (jCassam-Chenai et al.ll2004l) . In addition, we hnd that the peak position of Ne-Iike Fe 
L lines (0.7-0.85 keV) is located systematically inside that of the more ionized Fe L lines 
(0.9-1.2keV). This fact suggests that in each remnant, the Fe is indeed primarily SN ejecta 
heated by a reverse shock propagating toward the center. These radial prohles are depro- 
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jected by using the “onion-peeling” method, which estimates the emi ssion at each radiu s 
repeated removal of the contribution from concentric outer shells (e.g., iFabian et al.lIlQSOh . 
The deprojected prohles are displayed in the lower panels in Fig. 01 where we do not give 
error bars as they strongly depend on the geometry of the ejecta, but the values should be 
< ±10% if the assumed spherical symmetric distribution of the ejecta is correct. 

These deprojected prohles are used to estimate emitting volumes of the individual vnei 
components. Since the prohles of Kepler’s SNR look similar to those of SNR 0509-67.5, 
we assume that the emitting geometries are the same. The reverse shock radii, Rrs, are 
taken to be the sharp rises of the Fe K (and Ne-like Fe L), which are 0.7 and 0.65 times the 
forward shock radii, Rfs, for Kepler’s SNR and Tycho’s SNR, respectively, where we dehne 
Rps from eye inspections to be 107", 256", and 14.8" for Kepler’s SNR, Tycho’s SNR, and 
SNR 0509-67.5, respectively. We note that the value of Rrr found in this w ay for Tycho’s 
SNR is between those from a previous Chandra study ( Warren et ah 12005) an d a recent 


measurement from Fe K/3’s radial prohles with Suzaku flYamaguchi et ah 


2014). The hot. 


low-ionization, Fe-rich plasma (ejecta 3 in Table |2]) resides in a layer from Rrs to 0.85 (0.75) 
RpS) he., the inner radius of the half maximum of Si K’s deprojected prohles for Kepler’s SNR 
(Tycho’s SNR). Other IME-rich components, i.e., ejecta2 and ejectaS in Table |2]), occupy 
a shell outside the inner Fe-rich layers. The outer extent of the IME-rich components for 
Kepler’s SNR (Tycho’s SNR) is taken to be 0.97 (0.93) Rps, which corresponds to the position 
where we see a sharp rise of the deprojected Si K prohle for Kepler’s SNR (Tycho’s SNR). 
We note that the value for Tycho’s SNR repre sents the contact discontinuity suggested from 
a previous Chandra study flWarren et al.ll2005l) . We allocate volume hlling factors for the two 
IME-rich vnei components, such that their electron pressures are equalized: (/eji : /ej2) = 
(0.05 : 0.95), (0.2 : 0.8), and (0.015 : 0.985) for Kepler’s SNR, Tycho’s SNR, and SNR 0509- 
67.5, respectively. We then compute densities and masses for individual elements in each 
ejecta component. 


We note that the surface brightness of Kepler’s SNR is quite asymmetric; the northern 
rim is much brighter than the southern rim. This could lead to biased mass estimates, if we 
use only one side (either north or south) of the rim in estimating emitting volumes. There¬ 
fore, we examine radial prohles of the southern portion in addition to those of the northern 
portion. Ending that the reverse shock position as well as the IGE and IME layers are close 
to those of Tycho’s SNR rather than the northern portion of Kepler’s SNR. However, this 
modihed emitting volumes result in only ~10% variations of ejecta masses toward increas¬ 
ing/decreasing direction for IGE/IME. Such slight differences will not change our conclusion 
given below. 


The estimates of integrated masses for the three ejecta components are summarized 
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in Table |3l We find that IGE masses for SNR 0509-67.5 and Kepler’s SNR are close to 
each other, but are higher than that for Tycho’s SNR. This result quantitatively confirms 
our expectation (from the spectral comparison in Figs. [U and [3]) that Kepler’s SN was an 
overluminous (91T-like) la event. This is a lso consistent with th e result drawn by hydro- 
dynamical modeling of the X-ray spectrum fjPatnaude et al.ll2012l) . A more convincing plot 
can be found in Fig. [6l i n which we plot our data points on the IGF-IMF mass diagram of 
23 extragalactic SNe la (iMazzali et al.l 120071) as in Fig. [6l where several 91T-hke SNe are 
clustered in the upper left corner, a few subluminous (or 91bg-hke) objects are distributed in 
the lower right corner, and many normal la objects are located in between them. We can see 
that Kepler’s SNR is well within the overluminous (91T-like) category, while Tycho’s SNR is 
among the normal SNe la. In this plot, SNR 0509-67.5 can be found in a marginal region be¬ 
tween overluminous and normal la events. Threfore, this plot shows that the amount of IGE 
of Kepler’s SNR is even larger than that of the confirmed overluminous la SNR 0509-67.5, 
strongly supporting the idea that Kepler’s SN was an overluminous la event. 

It should be noted in Table [3]that (1) we provide fairly conservative uncertainties which 
are much larger than those from the statistical uncer tainties in Table El (2) we list unshocked 
IGE masses, based on SN la nucleosynthesis models (INomoto et al.lll984 iMaeda et al.ll2010[) 
combined with the reverse-shock positi ons estimated above ; and (3) we assume that the total 
IME -|- IGE mass becomes ~1.O5M0 flMazzali et al.l 120071) by adjusting the distances to the 
remnants. The details are described below. 


To evaluate the errors, we consider a wider range of the electron temperature for the hot 
Fe-rich component than the statistical uncertainties, given possible systematic uncertainties 
on Fe L emissivities. In fact, even if we fix fcTe =2keV and kT^ =15keV for this component, 
we derive similar quality fits; reduced-y^ values for kT^ =2keV and 15keV are, respectively, 
1.54 and 1.49 for Kepler’s SNR, 3.55 and 3.02 for Tycho’s SNR, and 1.42 (for both 2keV 
and 15keV) for SNR 0509-67.5. Thus, we consider these extreme temperature cases to be 
viable, and they provide the major source of the large uncertainties in Table [3l 


The upper two rows of Figure [5] show the masses of shocked IME and IGE as a func¬ 
tion of distance for three cases: the best fits in the red solid curves; kT^ = 2keV in the 
black dashed curves; and kT^ = 15keV in the black dotted curves. The third row shows 
unshocked eje cta masses for two S N la abundance distributions: the classical W7 model (the 


upper curve: INomoto et al.l 119841) and a two-dimensional (azimuthally -averaged') delayed - 
detonation model, i.e., off-center ignition 0-DDT model (the lower curve: iMaeda et al.ll2010l) . 
While both of the two models are not responsible for overluminous la SNe but for normal la 
SNe, the amount and composition of the unshocked ejecta would not be different by much, 
since the abundance distributions of the two models examined here look (more) similar to 
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those of overluminous cases rather than subluminous cases bv ISeitenzahl et al.l (120131) . To es¬ 
timate the cold ejecta masses, we integrate the ejecta masses from Vej = 0 np to = i?rs/age 
(the free expansion velocity at the reverse shock). We hnd that the composition of the un¬ 
shocked ejecta is almost pure IGE for both models, hence we name the y-axis label as “cold 
IGE.” The bottom panels show the total masses (hot IME -|- hot IGE -|- cold IGE). 

In principle, the ejecta masses can take any values (even unrealistic ones), by changing 
the temperature and/or distances (see Fig. E]). Therefore, we introduce one constraint that 
the total masses of IGE plus IME are within 0.96-1.14 M^t^, according to a sy stematic analysis 
of 23 samples of well-observed extragalactic SNe la flMazzali et al.N2007() . In the bottom 
panels, this mass range is given as horizontal dotted lines. Limiting the total masses in this 
range, we can constrain distances to the remnants for each temperature case. The distances 
for the best-ht cases (in between the two solid red lines) are ~4.2 kpc, ~2.55 kpc, and ~46 kpc 
for Kepler’s SNR, Tycho’s SNR, and SNR 0509-67.5, res pectively. While the values for the 


two Galactic SNRs agree with some of previous estimates (ISankrit et al.ll2005l: iKatsuda et ah 


20081: IVinkll2008l: iGhevalier et al.lll980 : Albinson et a 


smaller than the canonical values (IKerzendorf et ah 


2014 


19861) . all of the distances are slightly 


Havato et al.l I 2 OIOI: lAlvesI 12004 


Macri et al.1120061) . In particular, the distance to SNR 0509-67.5 in the LMG is re latively well 
determined. Thus, if we set the distance to the current best estimate of 48.1 kpc flMacri et ah 


2006 ). we obtain the total IME-I-IGE mass to be ~1.2M0 (s ee, Fig. |5l). This is j ust slightly 
larger than the upper end (1.14 Mq) expected for SNe la ( Mazzali et ah 20071) . and such 
a deviation would be still acceptable. On the other hand, there is also a possibility that 
SNR 0509-67.5 whose absorption within the LMG is negligible is indeed located at ~2kpc 
front side with respect to the main body of the LMG. Since accurate measurements of the 
total masses and distances, which require significant additional analysis efforts, are beyond 
the scope of this paper, we leave this issue for future work base d on high-resolution X - 
ray spectra which will be available with the upcoming ASTRO-H flTakahashi et al.l 120141) . 
Keeping the total IME-I-IGE masses within 0.96-1.14 Mq, we calculate allowed masses of the 
IME and IGE for the three temperature cases {kT^ =2keV, 15keV, and the best-ht). The 
resultant nominal values along with the uncertainties are given in Table |3l 


4. The Circumstellar Medium in Kepler’s SNR 


We perform spatially-resolved spectral analyses to reveal details of the GSM, using 
the XMM-Newton RGS as well as the Chandra AG IS. Our particu l ar int erest is to search 
for X-ray emission from the “diffuse” GSM found by IWilliams et al.l (120121) using Spitzer IR 
observations in the form of a lower temperature dust component, distinct from the hotter dust 
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component associated with dense/knotty CSM. To this end, we extract an RGS spectrum 
from the region within the solid lines in Fig. [TJ where no CSM knots are present. Figure [8] 
displays the RGS spectrum together with the best-£t model — the same model as we applied 
to the integrated spectra in Fig. [2J We hnd a line feature at ~0.5keV, which can be modeled 
by N Lya arising from the swept-up component. The best-fit parameters are listed in Table [21 
they show the presence of an enhanced N abundance in the swept-up component, similar 
to that found for the entire remnant (Figure [2] and Table [2]). This is evidence that this 
component originates from the CN O-processed CSM, rev ealing that the diffuse CSM can be 
seen in X-rays as well as in the IR fjWilliams et al.ll2012[) . 


The diffuse CSM can be reasonably interpreted as material in steady outflow from the 
progenitor system, which rehlled a cavity created by the major eruption a few 10^ yrs ag o 
suggested by optical proper motions of the CSM knots feandiera fc van den Bergh 1991 1. 
The density and the mass of the diffuse CSM are calculated to be 2.8 cm“^ and 0.08 Mq, 
respectively, at a distance of 4.2 kpc, based on the assumption that the emitting region is a 
polar-cap shell of 0.97-1 Rps with a filling factor of un ity. This density is consistent with 
that inferred from IR observations ([Williams et al.ll2012[l . The total mass integrated over the 
entire remnant is ~0.3 Mq, given that the region we examined covers only ~30% of the entire 
remnant. It should be noted that the total mass might increase by a factor of a few if we take 
into account the north-south asymmetry seen in both the IR and X-rays. Assuming a wind 
speed of =10kms“^, we estimate a mass-loss rate of 1.5xl0“^(uw/10kms“^) Mq yr“^, 
is needed to provide the estimated mass of 0.3 Mq within a radius of 2.1 pc at a distance 
of 4.2 kpc. Such a large mass-loss rate suggests that the progenitor is an AGB star, which 


is consistent with re s ults from recent hydro dynamical simulations (IVelazquez et al.[ [2006 
Chiotellis et ah 2012 : Toledo-Rov et ah 2014). 


Next, we turn to the bright knotty features which are positionally coincident with op¬ 
tically bright knots. We examine five such features (Knots 1-5) as shown in Fig. [71 By 
subtracting local BGs (dashed regions next to the source regions), we extract clean CSM 
spectra. The spectra are fitted by a single-component vpshock model. In the fitting proce¬ 
dure, the kTg, Hgt, and normalization are treated as free parameters. Abundances of O, Ne, 
and Mg are allowed to vary freely for the brightest Knots 1 and 2, while we fix all metal 
abundances to the solar values for the other relatively-dim knots. The spectra with the 
best-fit models are shown in Fig. [9l The best-fit parameters are summarized in Table [H The 
temperatures and ionization timescales are roughly consistent with those d erived from the 


integ rated RGS spectrum (Table [2]) and previous Chandra measurements ([Reynolds et ah 
20071 1. 


To derive densities and masses, we assume the X-ray-emitting plasma’s geometry to be 
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spherical for Knot-2 and a rugby-ball shape for the others. We also assume hlling factors 
to be unity for all the regions. The estimated densities and masses are listed i n Table HI 
The densities are an order of magnitude lower than the optical measurements flDennefeld 
1982[l . suggesting that X-rays arise from relatively tenuous regions surrounding (or close to) 
optical dense cores. The total mass of the five regions is I.IxIO^^Mq. Clearly, these regions 
contain quite limited amounts of the total CSM knots, since the optical CSM knots ar e 
distributed more widely than the selected regions (e.g., Bandiera fc van den BerghI Il99lh . 
Indeed, the normalization summed-up for the five regions (see. Table 0]) is only ~l/40 of 
the RGS-measured normalization integrated for the entire remnant (see. Table [2]). Part 
of the integrated normalization should be attributed to the diffuse CSM, whose fraction is 
estimated to be ~20%, based on the normalizations of the entire remnant and the southern 
part in Table [2J Thus, we infer the total mass of the knotty CSM in the entire remnant 
to be ~0.035Mq (=1.1x10“^Mq x 40 x 0.8). This would be slightly increased by adding 
optically emitting knots. Such a large amount of the CSM is consistent with those observed 
in SNe la-CSM. 

Unlike la-CSM SNe, the dense/knotty CSM in Kepler’s SNR is located far away from 
the progenitor star, as evidenced by three independent pie ces of information. First, the slow 
proper motions of the CSM knots measured in optical fiBandiera fc van den BerghI Il99lfl 
allow movements of only 5.3" (or 0.1 pc at a distance of 4.2 kpc) after the explosion. This is 
evidence that the knots were already ~2pc (i.e., the SNR radius) away from the progenitor 
star at the time of the SN explosion. Second, dividing the ionization timescale by the 
electron density in Table 01 we infer the shock-heating time to be ~50-200 yrs. This means 
that interactions between the CSM and the blastwave started at a very late phase, i.e., 
>200 yrs after the explosion. Third, the light curve of Kepler’s SN does not show any hints 
of early-phase CSM interactions. Figure [10] shows two ligh t curves of Kepler’s SN obtained by 


European observers in black and Korean observers in red flBaadelll943l: [Clark fc Stephenson 


1977 : Schaefer 1996h . For comparison, we show ty pical light curves of sev eral kinds of SNe , 
i.e., la-CSM (SN 200 2ic: Wood-vasev et ah 2004 1. 91T-like (SN 1991T: Lira et ah 1998 1. 


91bg -like (SN 1991bg: iFilippenko et al.lll992l) . and Type Up (SN 1999em: lElmhamdi et al. 
20031) . All of the light curves are normalized such that their peak magnitude is zero, in 


order to highlight the difference of their shapes. We hnd that the light curves of Kepler’s SN 
(both European and Korean versions) are clearly distinct from that of SN 2002ic (la-CSM), 
showing that there was no massive CSM in close proximity to the progenitor (<0.01 pc). 
This assertion would be fairly robust, given that both of the two independent historical light 
curves, generally agreeing with each other within 1 mag, decline quite fast (Amag~5 at 
day 100 from the maximum light), which is in stark contrast to the flat light curves seen 
for SNe la-CSM. 
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5. Discussion 


Our work strongly suggests that Kepler’s SN was likely to be a 91T-like event, based on 
the ejecta abundances — the dominance of IGE. Since detailed light curves of Kepler’s SN 
were taken by both European and Korean observers in A.D.1604-1605, it is interesting to ask 
if these light curves are consistent with those of overluminous SNe la obtained by modern 
telescopes. Unfortunately, we found it difficult to discriminate the overluminous and sublu- 
minous SNe la from the light-curve shapes. As shown in Fig. [101 fh® light curve measured by 
Johannes Kepler and some other European is best represented by that of SN 1991T, whereas 
the Korean version (red data in Fig. [TU]) is closer to that of SN 1991bg (subluminous). This 
discrepancy probably means that brightness uncertainties from naked-eye measurements are 
too large to discriminate among SNe la subtypes that follow th e light-curve d ecline rate 


versus peak luminosity correlation (i.e., the “Phillips relation”: jPhillipsi Il999h . Another 


important piece of information from the light curve is the peak brightness. Assuming the 
absolute maximum brightness of Kepler’s SN t o be -19.5 mag , a typical value for 91T-like 
events, and a visual extinction to be 3.27 mag f)Schaefeiill996l) . we estimate the distance to 
Kepler’s SNR to be either 6.3 kpc or 4.5 kpc from the apparent maximum brightness (-2.25 
or -2.95) reported in Europe and Korea, respectively. Therefore, the apparent maximum 
brightness from the Korean observations matches our result better. We should, however, 
note that the distance of 6.3 kpc from European maximum brightness is more consistent 
with recently preferred value of ~6kpc. Given these considerations, the only solid conclu¬ 
sion from the light curve studies is that Kepler’s SN was different from la-GSM and Type 
Up, but was close to Phillips-relation SNe la. 


In addition, we confirmed the presence of dense X-ray emitting GSM in Kepler’s SNR, 
and estimated its total mass in the entire remnant to be ~0.035 Mq. This large mass reminds 
us of the massive GSM associated with SNe la-GSM. However, as we discussed above, the 
light curve is inconsistent with a la-GSM event for Kepler’s SN, suggesting the different 
distances between the GSM and the progenitor star; a few pc away from the progenitor for 
Kepler’s SN is much farther than those seen in la-GSM events (<0.01 pc). Thus, we argue 
that Kepler’s SN was likely to be a 91T-like event, and started to interact with the dense 
GSM >200 yrs after the SN explosion. 


In this way, Kepler’s SNR has provided us with observational evidence for a recently 
proposed connection between 91T-like events and la-GSM objects: i.e., the two la subclasses 
are actually the same phenomeno n with seeming differe nces due to different timing of the 
interaction with a massive GSM (ILeloudas et ah! 1201311 . This speculation was originally 
inspired by optical observations of SN 2002ic and PTFllkx, whose optical spectra initially 
resembled a 91T-hke template, but started to develop an Ha line after 22 days and 59 days 
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after explosions, respectively, evolving to more like a Type Iln template. Since the timing of 
the CSM interaction is not nnique, we can easily gness a possibility that the timing was too 
late to show signatnres of the CSM, as the SNe fade after the explosions. While snch SNe 
would be classihed as 91T-like events, if we could monitor them for several hundred years, 
we should be able to detect a signature of a massive CSM. We believe that Kepler’s SN must 
be such an example, and speculate that another overluminous la remnant, SNR 0509-67.5, 
might show strong interaction with a massive CSM in t he future. It shoul d be noted that, 
whereas la-CSM objects occupy only OT -1% of SNe la flDildav et al.l 120121) . the fraction of 
91T-like events is ~9% flLi et al.ll2011al) . Therefore, our result implies that roughly 10% of 
SNe la are associated with a massive CSM. 

Based on the observational clues obtained so far, we discuss possible progenitor systems 
for Kepler’s SNR. First of all, the presence of the CSM i tself favors the SD scenario. While 
some recent DD models pr edict the existence of CSMs fjShen et al.l 120131: lRa.skin fc Kasen 
20131: ITanikawa et al.ll2015l) . they originate from a H-rich envelope of a He WD or tidal 
tails from a C/0 WD, hence they would not be CNO-processed N-rich material. This is 
in conflict with the overabundance of N seen in the CSM of Kepler’s SNR. Moreover, the 
small amount of O ejecta, <0.069 M^n as sho wn in Table El is not expected for DD scenarios 
fe.g.. IPakmor et al.ll201ll: iRopke et al.ll2012lj. but is mo r e consistent with delaye d-detonation 


Chandrasekhar-mass models (e.g., iMaeda et al.l l2010l: ISeitenzahl et al.l 120131) . Therefore, 


there is no observational support for DD channels as the origin of Kepler’s SN, except for 


the absence of a candidate surviving companion star flKerzendorf et al 


201411 . 


Outflows from SD progenitors can originate from either the exploding star or the donor 
star. For Kepler’s SN, the former case is less likely than the latter. This is because the 
exploding star, which must have experienced a major mass loss ~ 3 x 10"^yrs before the 
explosion, cannot have had sufficient time to evolve to a WD and accrete material from 
its donor star to near the Chandrasekhar mass (cf., accretion timescale is of the order of 
lO^yrs). On the other hand, the exploding star is not necessarily a WD — the so-called 
Type 1.5 model predicts that a degenerate C/0 core of a massive ABC star wi th very 


low metallicity grows to near the Chandras ekhar mass before it loses its envelope flArnett 


19691: llben fc Renzinil Il983l: ILau et al.l 1200811 . However, such a SN occurs inside a massive 
H envelope, so that the SN light curve should look like Type Hp SNe as shown in Fig. [TU] 
(K.Maeda et al. in preparation). Looking at Fig. [TOl this is evidently not the case for Kepler’s 
SN. Moreover, the fact that the metallicity of Kepler’s CSM is nearly solar (Table H]) also 
conflicts with the Type 1.5 model. 

Turning to outflows from a companion/donor star, the high mass-loss rate we estimated, 
M ~ 1.5 X 10“®(nw/10kms“^) MQyr“^, suggests an AGB companion. One possible progeni- 
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tor system hosting an AGB companion is a recurrent nova flHachisu &: Katoll200ll) . However, 
there are two major observational facts contrary to the recurrent novae scenario. First, Ke¬ 
pler’s SN experienced a major mass eruption ~ 3 x lO'^yrs before the SN explosion. Such 
a long nova timescale is not expected for recurrent nova systems with near Chandrasekhar- 
mass WDs, f or which typical recur rence times are 10-100 yrs. Second, according to numerical 
simulations (IMarietta et al.l 120001) . the impact of SN ejecta removes most of the envelope of 
a companion star if it is a giant like an AGB star. As a result, the surviving companion will 
evolve away from the red giant branch, along a track of constant luminosity ( r'U 10^ Lq) for at 
least 10® yrs. Such a bright star, if present, should be easily d etected by optical obser vations, 
but is not detected at least near the center of the remnant flKerzendorf et al.l 120141) . Thus, 
a recurrent nova is not likely the progenitor of Kepler’s SN, either. 

There are other models that satisfy observational properties of Kepler’s SNR, i.e., both 
high mass-loss rate and no bright surviving cornpanion . In efforts to explain the progenitor 
system of SNe la-CSM, iHan fc Podsiadlowski (1200611 investigated the classical supersoft 
channel that consists of a C/0 WD and a main sequence star. They found that if the donor 
star is a relatively massive (^SMq) main sequence star and the system experiences delayed 
dynamical instabilities, a large amount of mass can be lost from the system in the last 
few 10"^ yrs before the explosion. As a result, the donor star loses most of the mass before 
the explosion, and its luminosity becomes as low as the sun. This model can explain both 
the massive GSM and no detectable surviving companion in Kepler’s SNR. Alternatively, 
Soker et al.l (1201311 argued that a WD and a hot core of a massive AGB star in a common 
envelope will violently and promptly merge if the WD is denser than the core. In this 
case, there is no surviving companion left behind the SN while massive a GSM could be 
present near the progenitor system, making this model viable as well. However, we have 
not yet found bservational evidence supporting these models, leaving us the mystery of the 
progenitor system. Continuous searches for a surviving companion and/or signatures of a 
companion star in Kepler’s SNR would be of great importance for revealing the progenitor 
of Kepler’s SNR and SNe la in general. 


6. Conclusions 

We have analyzed X-ray spectra and images of three young (~400yrs old), ejecta- 
dominated SNRs, i.e., Kepler’s SNR, Tycho’s SNR, and SNR 0509-67.5 in the LMC. A simple 
spectral comparison among the three remnants suggests that Kepler’s SNR and SNR 0509- 
67.5 produced substantially more IGE than Tycho’s SNR. In fact, we quantitatively revealed 
that the IGE masses for both Kepler’s SNR and SNR 0509-67.5 are similar to those of 
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overluminous SNe la, while that for Tycho’s SNR is com parable to that of ii ormal SNe la, 
assuming the total mass of IME and IGE to be ~1.05 Mq flMazzali et al.ll2007() and inferring 
the unshocked ejecta masses using the SN la’s nucleosynthetic models together with the 
reverse shock positions estimated from Chandra images. This inference is consistent with 
the fact that SNR 0509-67.5 and Tycho’s SNR are confirmed remnants of overluminous and 
normal SNe la explosions, respectively. Therefore, we conclude that Kepler’s SN was likely 
an overluminous (91T-hke) event. Our method is quite simple, and thus could be a powerful 
tool for discriminating subclasses of Type la SNRs. 

The XMM-Newton's RGS spectrum shows a strong N Lya line for Kepler’s SNR. We 
measured the N/H abundance to be ~4 times the solar value, suggesting the presence of a 
GNO-processed GSM as was previously suggested by optical spectroscopy flBlair et al.lll99lf) . 
We revealed that the GSM is comprised of diffuse gas and dense knots. The GSM masses are 
estimated to be ~0.3Mq and ~O.O35M0 for the diffuse gas and dense knots, respectively. 
The mass of the diffuse GSM allows us to infer a mass-loss rate of the progenitor system 
to be ~ 1.5 X 10“®(nw/10 kms“^) Mq yr“^, consistent with those of AGB stars. As for the 
knots, their slow optical proper motions and ionization timescales measured by using X-ray 
spectra suggest that they were located far away from the progenitor system at the time of SN 
explosion, and were recently shock heated. We argue that Kepler’s SN was likely a 91T-like 
event that started to interact with massive GSMs a few hundred years after the explosion. 
This supports the possible link between 91T-like SNe and SNe la-GSM, implying that ~10% 
of SNe la are associated with massive GSM. All of these SNe would be related to a single 
degenerate channel. 
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Energy (keV) 


Fig. 1.— Suzaku XIS spectra of Kepler’s SNR in red, Tycho’s SNR in black, and SNR 0509- 
67.5 in blue. The intensities are normalized such that the peak values of Si K-shell lines are 
equalized. The error bars are omitted for clarity. The Fe L and Fe K lines for Kepler’s SNR 
and SNR 0509-67.5 are clearly stronger than those of Tycho’s SNR. It should be noted that 
these spectra are raw data, and are not corrected for the interstellar absorption. 
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Fig. 2.— Spatially integrated X-ray spectra for Kepler’s SNR (upper panel), Tycho’s SNR 
(middle panel), and SNR 0509-67.5 (lower panel). We use the XMM-Newton RGS below 
2keV and the Suzaku XIS above 2keV for both Kepler’s SNR and SNR 0509-67.5, but only 
the Suzaku XIS for Tycho’s SNR. For Kepler’s SNR, we simultaneously £t the Chandra ACIS 
spectrum extracted from the GSM knots (Knots 1 and 2 in Fig. [7]). The lower panels show 
residuals. Individual components, i.e., swept-up, three ejecta, and power-law, are separately 
shown in blue, red-|-orange-l-light blue, and gray, respectively. 
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Fig. 3.— Comparison of unabsorbed ejecta emission models for Kepler’s SNR in red, 
SNR 0509-67.5 in blue, and Tycho’s SNR in black. The models are convolved with the 
XIS response function, and the intensities are normalized such that the peaks of Si K-shell 
lines are equalized. The lower panel shows spectral ratios with respect to Tycho’s SNR. 
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Fig. 4.— Top panels: Chandra ACIS images of the three SNRs analyzed in this paper. 
We extract radial prohles from the white pie-shaped areas covering relatively smooth, limb- 
brightened regions. Middle panels: X-ray radial prohles from the pie-shaped regions in top 
panels. The prohles are generated in four diherent energy bands: 0.7-0.85keV (black), 0.9- 
1.2keV (red), 1.65-2.0keV (green), and 6.1-6.8keV (blue). Bottom panels: Deprojected 
radial prohles corresponding to the prohles in the middle panel. 
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Fig. 5.— Ejecta masses as a function of distances to Kepler’s SNR (left column), Tycho’s 
SNR (central column), and SNR 0509-67.5 (right column). The top-row, second-row, third- 
row, and the bottom-row displays masses for shocked IME (Si-|-S-|-Ar-|-Ca), shocked IGE 
(Fe-|-Ni), unshocked IGE, and their sum, respectively. The dashed and dotted curves are 
obtained by hxing the electron temperature of the Fe-K emitting plasma at 2 keV and 15 keV, 
respectively, whereas solid red curves are responsible for the best-£t results. The cold IG E 


mass es are based on two nucleosynt hetic models, i.e., W7 (the upper curve: iNomoto et ah 


1984J ) and 0-DDT (the lower curve: iMaeda et al.ll2010l ). These two cases result in the pairs 


of curves in the bottom panel. In the bottom panels, the horizontal lines indicate the range 
of O.96-1.14M0 Mazzali et ah (2007). 
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Fig. 6.— IGE masses vs. IMF masses for 23 extragalactic SNe la (crosses: iMazzali et al. 


20071) . The data include a variety of subclasses of SNe la, i.e., overluminous (e.g., SN 1991T), 


normal (e.g., SN 1994D), and subluminous (e.g., SN 1991bg). Our data for Kepler’s SNR, 
Tycho’s SNR, and SNR 0509-67.5 are shown as a red triangle, a black box, and a blue 
diamond, respectively. The errors quoted are taken from Table El 
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Fig. 7.— Spectral extraction regions for spatially-resolved CSM studies in Kepler’s SNR. For 
the diffuse CSM region, we extract the XMM-Newton^s RGS and Chandra^s ACTS spectra 
from the region within two white lines. For the CSM knots, we extract Chandra’s ACIS 
spectra from a solid circle and four solid ellipses. Dashed regions are where we extract local 
BGs. 
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Fig. 8.— The combined RGS (below 2keV) and the Chandra's ACIS (above 2keV) spectrum 
extracted from the southern portion of Kepler’s SNR as shown in Fig. [71 The spectrum is 
htted with the same model as we applied to the integrated spectra in Fig. O The lower 
panels show residuals. 















32 



Energy (keV) Energy (keV) 



0.5 1 2 


Energy (keV) 


Fig. 9.— Local-BG subtracted Chandra's ACIS spectra for CSM knots as shown in Fig. [71 
The spectra are htted with a single-component vpshock model. The lower panels show 
residuals. 
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Days from maximum light 


Fig. 10.— Light curves of Kepler’s SN and other modern SNe. The triangular data points 
in black and red tied by solid lines are responsible for Kepler’s SN, obtained by Europea n 
and Korean observers, respectively. The errors are assumed to be ±0.25 mag (jBaadelll943h . 
Other hlled-circle data tied by dotted lines are responsible for recent SNe including la-CSM 
(SN 2002ic), 91T-like (SN 1991T), 91bg-like (SN 1991bg), and Type Up (SN 1999em). 











Table 1. X-ray data used in this paper 


Target 

Instrument 

Obs. ID 

Obs. Date 

Effective exposure (ks) 

Kepler’s SNR 

XMM-Newton^s RGS 

0084100101 

2001-03-10 

32.6 

Kepler’s SNR 

Chandra^s ACIS 

6715 

2006-08-03 

159.1 

Kepler’s SNR 

Suzaku's XIS 

505092040 

2011-02-28 

146.2 

Tycho’s SNR 

Suzaku's XIS 

500024010 

2006-06-27 

101.1 

SNR 0509-67.5 

XMM-Newton's RGS 

0111130201 

2000-07-04 

36.0 

SNR 0509-67.5 

Suzaku's XIS 

508072010 

2013-04-11 

175.9 


Note. — The Pis of the data are A. Decourchelle (Kepler - XMM), S. Reynolds (Kepler 
- Chandra), S. Park (Kepler - Suzaku), Suzaku science working group (Tycho - Suzaku), 
M. Watson (SNR 0509-67.5 - XMM), and H. Yamaguchi (SNR 0509-67.5 - Suzaku). 
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Table 2. 


Spectral-fit parameters for the integrated X-ray spectra 



Parameter 


Kepler’s SNR 

Kepler’s SNR (south) 

Tycho’s SNR 

SNR 0509-67.5 

CSM component 







kTe (keV) 


1.06±0.03 

0.4+0.20 

0.41+0.01 

1 76+6-6® 

-1- '0-0.48 


log(net/cm~^ sec) 


lO-sifo.S? 

10.651“;““ 

xu.oy_Q Q 5 

10 22+'^ '^® 


Abundance‘s (solar) N 


Q Q-1 +0.24 
'^■'^-'^-0.25 

6.08t“;" 

i*> 

0.135“ 


Redshift (10“^) 


1.39±0.05 

=Ejectal 

= Ejectal 

=Ejectal 


Line broadening (E/lkeV eV) 


9 q-i +U.CJD 

=Ejectal 

= Ejectal 

=Ejectal 


f neniidV/4:7Td‘^ (lO^^crn”®) 


I93.i5;^;“| 

9.37;+p 

1788.321!+" 

1 55+6-14 

1.00 Q 43 

Ejecta components 






1 

kTe (keV) 


0.37±0.01 

0.44l“o:S^ 

o.7ol!;" 

1 09+°-^^ 
-^•6’^-0.24 


log(net/cm~^ sec) 


10.52+0.01 

1A oq+0.06 

iU.3y_o_o3 

10.79+0.01 

10 14+6-02 
XU.X‘i_Q 01 


Abundance (10'^ solar) C 


0.54*’ 

0.54*’ 

0.54'> 

0.54“ 


N 


o'’ 

o" 

o'’ 

o'’ 


O 


0 9K+0-02 

U -^0_Q Q3 

0.31+:m 

0.75+0.01 

0 22+'^'^ 
6-^^-0.15 


Ne 


0 67+0-0® 
0-D<_o,o7 

0.84+;!^ 

0(< 0.0003) 

-,=+0.18 

l-lo_o.i7 


Mg 


0 77+0-0^ 

' ' -0.06 

0.53+0.14 

0.34+0.01 

6-6'-o.07 


Si 


lo*- 

lO*- 

lo’’ 

io“ 


s 


15.53+:“ 

12 31 + 0-64 
iz.qi_o eg 

10.651!;" 

15.791!;!! 


Ar 


18.98t°;« 

11 34+1-00 
xx.o^_l_7g 

11.021!;" 

10.591!;!! 


Ca 


36 40'''^'^^ 

OD.‘iU_^ gQ 

23 74+®-7-^ 

'‘•5- '^-5.62 

25.831!;!! 

14 59+*^'® 
x-i.oy_2 46 


Fe 


13.65+;" 

6.63f+j? 

1.49+0.01 

175.561!!-““ 


f rieniidV/ATvd'^ (10^cm~^) 


961.42j:ll|‘)2 

230.29tl“+f 

9014.62l45^“g, 

0.i3l!;!! 


Redshift (10“^) 


-2.94+0.01 

-5.6oli;“t 

-4.57+0.01 

-0 89'^'^'^^ 
u.oy_Q 4 ^ 


Line broadening (E/lkeV eV) 


Q -I o + O. 14 
’’^■-'^'^’-0.21 

10 16+0-*^® 
xu.xo_o_45 

».»u_o 03 

11.901!;" 

2 

kTe (keV) 


9 nR'*”0-0i 

2.U«_q Q2 

-'--‘^D-O.OS 

0.96+0.01 

5.441!;" 


log(net/cm~^ sec) 


10.32+0.01 

=Ejectal 

10.90+0.01 

=Ejectal 


Abundance‘s (10^ solar) Fe 


3 co+0.04 
05.O8 _q 05 

Q Qo+0.15 

0(< 0.0002) 

0(< 0.03) 


f neUiidVfAird'^ (10®cm~®) 


S75.35t.tll 

131.51++ 

22.75+0.02 

0.013+0.001 


Redshift (10~^) 


=Ejectal 

=Ejectal 

= Ejectal 

=Ejectal 


Line broadening (E/lkeV eV) 


=Ejectal 

=Ejectal 

= Ejectal 

=Ejectal 

3 

kTe (keV) 


2.59+0.01 

1.92+0.1 

9 34+6-63 

6 97+^'*^ 

o-y'_i,08 


log(net/cm~^ sec) 


9.21+“ 

q 02+0-04 
y-U-i_Q_Q3 

9.80+0.01 

y.ou_Q Q 2 


Abundance*^ (10“^ solar) Ar 


0 (< 0.19) 

0.45+0-^® 

0 (< 2.7) 

37 61+^'^'^ 


Ca 


1 43+0-31 

1 65+°-®9 

l.DO_o 58 

64.821!;“! 

47.281!“;!! 


Fe 


lo*- 

lo*- 

10“ 

io“ 


Redshift (10“^) 


-5 73'*''^'^^ 

=^-2O_i_09 

-5 24+6-66 

o.32l?:S 


Line broadening (E/lkeV eV) 


12.10+0.27 

11.151+1 

10 76+6-23 

XU. 'D_o,24 

16.811!;" 


f ne'niidVlATTd'^ (10®cm~®) 


5808.471^1;" 

1748.8011“+! 

1.08+0.01 

U.UUD_Q QQ^ 

Power-law component 







r 



9 ci:: + 0.02 

.^.OO_0_0l 

2.56+0.01 

q q7 + 0.16 
'^’•'^’'-0.2 


Norm (phkeV~^ cm~^ s~^ at 1 keV) 


98.93t5;51 

24 49+0-43 

.i-i.-iy 0 30 

714.681!;!! 

qb:+0.37 

1.»9 q 4g 

Additional lines 







Fe L-l-O K Center (keV) 


0.708t“;““ 

— 

0.741+0.001 

0.724+0.002 


Norm (10~^ phkeV~^cm“ 

^s-l) 

es.i+i.o 

— 

986.88+19.1 

2.38+0.3 


Fe L-|-Ne K Center (keV) 


— 

1.201+0.011 

• 0.003 

— 


Norm (10“^ phkeV“^cm~ 

2s-l) 

— 

1.96+! 

33.89+0.6 

— 


Fe L-|-Ne K Center (keV) 


1 997+0.001 
-^•^^'-0.003 

— 

1.254+0.001 

1.244+0.001 


Norm (lO”"^ phkeV“^ cm~ 

2s-l) 

22.75+ + 

— 

42.40+0.5 

0.19+0.1 


Cr K Center (keV) 


5 514+0-035 
O.014:_o 037 

5.4381!;"“ 

5 474+6.016 

^•^'^-0.019 

5.4361!;!!! 


Norm (10“^ phkeV”^ cm~ 

2b-1) 

85.18+^?'^ 

17.96+10.2 

197.98 

7.05+3.2 


Mn K Center (keV) 


5 076+0-04 

0.y<'D_o_o42 

c 000+0.266 

5 935+6-07 

O.yqo_o_o71 

5 cqq+0.068 

O.ey3-0 071 


Norm (10~^ phkeV”^ cm~ 

2s-l) 

60.18+16.3 

14 16+^0.8 
x-i.xo i4 2 

41.07+32 

3.56+3.4 


/d.o.f. 


8289.8 / 5724 

1744.2 / 1305 

5659.0 / 1981 

1031.6 / 783 


Note. — Nil values are fixed to 6.4x10^^ cm~^, 1x10^^ cm~^, and 0.6x10^^ cm~^ fo r Kepler’s SNR, Tycho’s SNR, and SNR 0509- 
67.5, res pectively. Other elements are fixed to ei ther the solar values IWilms et al.ll2000l^ or typical ISM values of the Large Magellanic 
Clouds iKorn et aljf2002| ; [Russell &: DopitaHl99^ for the two Galactic SNRs or SNR 0509-67.5, respectively. ^Fixed values. ‘^Other 
elemental abundances are lined to those in the Ejectal component. ‘^Other elemental abundances are fixed to zero. 
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Table 3. Ejecta masses (Mq) 


Element 

Kepler’s SNR 

Tycho’s SNR 

SNR 0509-67.5 

0 

0 01 

U.U10_g QQg 

0 31Q+0-012 
u.oiy_0 159 

U.UO^_0 018 

Ne 

U.UU^_o 004 

0(< 0.001) 

O 

O 

CO 

CO 

1 + 

o o 

b b 

to o 

CO to 

Mg 

0.002(< 0.005) 


0.001(< 0.002) 

Si 

U.U^U_o 014 

0.28AZZ 

U.iUD_o 072 

S 

-0.016 

o.zzy_o 096 


Ar 

0 01 

U.UiZ_o 006 

U.U/Z_o 024 

omsZZi 

Ca 

U.UZi_o 012 

nil 0+0.012 

U.ilQ_o_o41 

U.UOy_o 072 

Ee (shocked) 

0 7in+0-094 

U. 1 iU_o 249 

0 IIA+O.SS? 

U.iO4_0 028 

n Q79+0.355 

U.Of Z_o 118 

Ni“ (shocked) 

0 110+0015 
U.iiU_o 038 

0 021+0-060 
U.UZi_o 005 

omsZZ 

IGE^ (unshocked) 

u.izy_o 087 

0 104+0-099 

u.ia4_o 120 

0.318l°;i^ 


Note. — “The Ni /Fe abundance ratio is assumed to be that expected 
in the W7 model ( Nomoto et ah 1984 '). ^Taken from the cold IGE 
panels in Fig. [5l 


Table 4. Spectral-fit parameters for the GSM knots in Kepler’s SNR 


Parameter 

Knot-1 

Knot-2 

Knot-3 

Knot-4 

Knot-5 

fcTe (keV) 

l.lOzhO.ll 

0 95+°'^ 

u.yo_o 08 

1 qo + 057 

0-68t°:J| 

i-8R“:r, 

log(net/cm~^ sec) 

iu.oy_Q Qg 

iu.yD_Q j^2 

10 66"*'^'^^ 
1U.DD_q J^2 

lO-aRo.l? 

10.367“;=^ 

Abundance (solar) N 

3.5“ 

3.5“ 

3.5“ 

3.5“ 

3.5“ 

O 

0 97+°-^® 
-0.13 

u.yo_o_og 

1 

1 

1 

Ne 

1 “11 'hb-23 
-‘-•'^-‘--0.22 

'4-o^_0.21 

1 

1 

1 

Mg 

-^•^'^-0.18 


1 

1 

1 

f rieniidV/A-Trd'^ (lO^^cm”^) 

T 07+0.17 

1.58+0.13 

U.oo_o_o5 

0 77+0-14 
‘ ‘ -0.12 

tJ-^O_0.02 

Tie (cm~^) 

23.4 

29.2 

6.4 

11.3 

6.9 

Mass (10“'* M^) 

2.8 

1.9 

3.1 

2.4 

1.3 

X^/d.o.f. 

181.8 / 134 

163.1 / 111 

130.1 / 128 

95.9 / 109 

95.6 / 101 

Note. — Nil values are fixed 

to 6.4x10^^ cm 

We assume 

Tie = 1.2tijj 

and the distance 

d = 4.2 kpc. 


Fixed values. 














